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APPROXIMATIONS OF THE WIENER SAUSAGE AND ITS 
CURVATURE MEASURES 

By Jan Rataj 1 , Evgeny Spodarev and Daniel Meschenmoser 

Charles University, Ulm University and Ulm University 

A parallel neighborhood of a path of a Brownian motion is some- 
times called the Wiener sausage. We consider almost sure approxi- 
mations of this random set by a sequence of random polyconvex sets 
and show that the convergence of the corresponding mean curvature 
measures holds under certain conditions in two and three dimensions. 
Based on these convergence results, the mean curvature measures of 
the Wiener sausage are calculated numerically by Monte Carlo simu- 
lations in two dimensions. The corresponding approximation formu- 
lae are given. 

1. Introduction. Although the path of the Brownian motion is a ran- 
dom fractal set, the boundary of its parallel neighborhood of radius r (also 
called the Wiener sausage) is a Lipschitz manifold and the closure of the 
complement to the parallel neighborhood has positive reach for almost all 
r > 0, at least in dimensions two and three; see [5]. Due to this fact, curva- 
ture measures of the Wiener sausage, including its volume, surface area and 
Euler-Poincare characteristic, can be introduced, as in [13]. Since the 1930s, 
much attention has been paid to the computation of the mean volume of 
the Wiener sausage; see [8, 16] and [1]. Recently, its mean surface area was 
obtained in [12]. The same result was re-established, using different tools, 
in [9]. Explicit formulae for other mean curvature measures of the Wiener 
sausage are still unknown. The first attempt to treat them was made by 
Last (in [9]), who reduced the computation of his support measures to two 
mean curvature functions. 

In this paper, we consider the curvature measures of the Wiener sausage 
in dimensions two and three. In the planar case, we show that the mean 
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Euler number is finite, by using a general property of the set of critical 
values of the distance function and the self-similarity of the Brownian mo- 
tion. Furthermore, we consider a.s. approximations of the Wiener sausage by 
polyconvex sets and show that the corresponding curvature measures con- 
verge almost surely and, in certain cases, also in mean. Finally, we attempt 
to study the mean curvatures of the Wiener sausages by simulation and 
find the corresponding approximation formulae. For that, we consider the 
approximating polyconvex sets whose curvature measures can be calculated 
numerically. Averaging over a large number of Monte Carlo realizations, the 
mean curvatures of these sets are estimated and the corresponding approx- 
imation formulae in two dimensions are found. Although the formulae for 
the mean area and boundary length are known explicitly (see the discussion 
in Section 7.1), they involve an integral of a functional of Bessel functions 
that has to be assessed numerically. Thus, apart from illustrating our con- 
vergence results stated in Sections 5 and 6, our approximation formulae for 
the mean area and boundary length are of practical interest. As for the mean 
Euler characteristic EVo(H r ) of the Wiener sausage, no analytical expression 
is known thus far. Thus, the approximation formula in Section 7.2 is novel 
and can be considered as a first step toward studying its properties. 

The paper is organized as follows. In Section 2, curvature measures and 
intrinsic volumes are introduced for polyconvex sets and the sets of positive 
reach. Section 3 contains some preliminaries on Wiener sausages. In Sec- 
tion 4, we show that the mean Euler number of the Wiener sausage is finite in 
IR 2 . The rest of the paper deals with approximations of the Wiener sausage. 
Namely, the convergence of parallel neighborhoods of compact sets in the 
Hausdorff metric is considered in Section 5. In Theorem 5.2, conditions are 
found under which the weak convergence of their curvature measures takes 
place. Polyconvex approximations of the Wiener sausage satisfying these 
conditions are described in Section 6. Finally, Section 7 contains the Monte 
Carlo simulation results for the mean curvature measures of the Wiener 
sausage. 

2. Intrinsic volumes and curvature measures. Let be the Lebesgue 
measure and H s the s-dimensional Hausdorff measure in K~ (see, e.g., [14]). 
Denote by /C the set of all convex bodies in Mr. Let 1Z be the family of all 
finite unions of sets from /C. It is sometimes called the convex ring and its 
elements are referred to as polyconvex sets. 

Let A © B be the pointwise sum of two sets, A and B, in M~. For the 
closed ball B = B r (o) of radius r > in M. d centered at the origin, the set 
A r = A@B r {o) is called the r-parallel neighborhood of A. The operation 
A i— > A r is known as dilation. 
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Let A be a nonempty compact subset of M. d . The distance function of A 
is defined by 

Aa : x m.in{\x - a\ : a e A} , xeR d . 

Given two nonempty compact subsets A, B of M. d , their Hausdorff distance 
is defined as 

dn (A,B) = max] max A# (a) , max (b)>. 

{ a&A b£B ' 

It is well known that du is a metric. For any x S R d , we denote by 

T,a(x) = {a E A : \x — a\ = A^(a;)} 

the set of all points in A which are nearest to x. The set is always 

nonempty, by the compactness of A. We say that x G M. d is a critical point 
of Aa if x lies in the closed convex hull of S^(x). A point x is called regular 
if it is not critical. A number r > is a critical value of Aa is there exists 
a critical point x of A^ with Aa(x) = r. We shall denote by C(A) C (0,oo) 
the set of all critical values of A^. 

The reach of a closed subset ^4 C M. d is given by 

reach A = sup{r > : Vx G M d , A^(x) < r =>• card Sa(^) = 1} 

(see [3]). As examples of sets with positive reach, consider convex closed sets 
(with infinite reach) or sets with compact and C 2 -smooth boundaries. 

For a compact set A C M. d with positive reach, the Steiner formula holds 
for sufficiently small radii. More exactly, we have 

d 

V d (A r ) =J2^yV d -i(A), < r < reach A, 

i=0 

where V d is the Lebesgue measure in M. d (volume), uj% is the volume of the 
unit i-ball and Vi is the ith intrinsic volume of A, i = 0, . . . ,d (see [3]). The 
functionals Vi are motion invariant and additive. In particular, Vq is the 
Euler-Poincare characteristic and Vd-\ is one half of the surface area. 

A local version of the Steiner formula holds for sets with positive reach. 
Let £a( x ) be the nearest point of A to x whenever A^x) < reach A Then, 
we have, for any bounded Borel subset F of M d , 

d 

V d ((A r \ A) n Q\F)) = J2^r l C d ^(A; F), < r < reach.4, 

where Ci(A; •) is the ith curvature measure of A; it is a signed Radon measure 
concentrated on dA for < i < d — 1 and C d -i(A;-) is the restriction of 
the (d — 1) -dimensional Hausdorff measure to dA, provided that A is d- 
dimensional. If dA is compact, then Ci(A; M. d ) = Vi(A), i = 0, . . . , d — 1. 
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The curvature measures Cj(-; •) admit an additive extension to the family 
of polyconvex sets; see [15]. The same holds, of course, for intrinsic volumes 
Vi(-), i = 0, ...,d. Curvature measures can further be extended to certain 
full-dimensional Lipschitz manifolds; see [13]. This applies, in particular, to 
the case when A is a full-dimensional Lipschitz manifold in M. d and K d \ A r 
has positive reach; see [13], Proposition 4. In this case, the curvature mea- 
sures of A r can be introduced through the identity 

Ci(A r ;F) = (-l)^- 1 ^^^; F) 

for i = 0, . . . , d — 1 and r > such that reach(M ci \ A r ) > 0. 

3. Mean curvatures of the Wiener sausage. Let (fi, J, P) be an arbitrary 
probability space. The Wiener process is a random process {VF(i) :t > 0} 
defined on (0, 3", P) with continuous paths and with the following properties: 

• V^(0) = almost surely; 

• W has independent increments; 

• W(t) — W{s) has Gaussian distribution with mean zero and variance t — s 
whenever < s < t. 

Given d independent Wiener processes W±, . . . ,Wd, the random function 
X(t) = (Wi(t), . . . , Wd(t)), t > 0, is the Brownian motion in M. d initiated at 
the origin o G W d ; see, for example, [2]. Denote by E = {X(t) : < t < 1} C M. d 
the image of the Brownian motion in M d running over < t < 1 . 

Definition 3.1. The set E r = S © B r (o), r > 0, is called a Wiener 
sausage; see, for example, [18], page 64. 

By [12], Corollary 4.4, the intrinsic volumes V^(H r ), i = 0, . . .,d— 1, are 
well defined for d < 3 almost surely. Moreover, the expectation EV^(H r ) is 
finite for i = d, d — 1; see [1] and [12]. It will be shown in the next section 
that EVb(H r ) < oo in the case d = 2. We conjecture that all mean intrinsic 
volumes of H r are also finite in three dimensions. 

4. Finiteness of the mean Euler number. In this section, we show that 
the mean Euler number EVo(H r ) exists and is finite for any r > in two 
dimensions. Two basic tools are used. The first is a general result about the 
smallness of the set C (A) of critical values of A C IR 2 , which is of a similar 
nature to the main result in [5]. Furthermore, the stochastic self-similarity 
of a Brownian motion is used to show a certain local homogeneity of the 
random set C(H). 

Lemma 4.1. Let 4CR 2 be compact with D = diam^4 and let < e < 
r be such that C{A) n [r, r + e) = 0. Then, the total number of bounded 
connected components of M? \ A r is less than or equal to -^^2=} . 
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Proof. Let L\,Li,--. be an enumeration of the bounded connected 
components of M 2 \^4 r and let £j be the inradius of Lj (i.e., e$ is the maximal 
number such that there exists a disc of radius e$ contained in Li). Let, for 
each i, Bi = B £i {si) be any fixed disc with this property. Then, Sj is a critical 
point of A^4 with critical value r + e$. 

Consider two different centers, Sj and s^. Since they belong to different 
connected components, the segment [si, Sj] must contain a point of A r . This 
point must belong to a disc of radius r which is disjoint with the interiors 
of both Bi and Bj . Simple geometric reasoning leads to 

I «i - sj\ > \J 2rEi + e 2 + <J 2rSj + e 2 . 

Let Zi = {z G R 2 : dist(z, Bi) < dist(z, B,-)Vj / be the "zone of influence" 
of Bi. We shall find a lower bound for the area of Zj. Let 6j be a boundary 
point of Bi belonging to A r ; thus, etj = 6j + ^(6j — Sj) G A Setting = 
6j + f:(si — bi), all three discs Bi, B r {ai) and -B £ (s^) have a common tangent 
at bi. Now, let i| and t 2 be two tangent points of dB r (ai) on lines passing 
through s[. It follows from the previous observations that the quadrangle 
aitjs'it 2 lies in Zi and its area is r\j2re + e 2 > r\j2re. Consequently, the 
total number of zones Z; must not exceed 7r (- p + 2r ) □ 

Lemma 4.2. iei n G N, £l, . . . ,£ n be positive numbers and So, ... , S n 
be segments of length at least s > contained in the ball Br(o), R> s. 
Assume that for allO <i < j < n, dist(Si, Sj) > \J s /2yjEi+i + • • • + Sj. Then, 
^+--- + ^< 8s~ 3 / 2 irR 2 . 

Proof. To prove the lemma, we shall use the following fact which can 
be shown by induction. Whenever n G N, 6±, . . . , 6 n > and Pq, . . . ,P n are 
points on a line such that 

dist({P , • • • , Pi-l}, {Pi, ...,P n })>8i, i = 1, ... ,71, 

then 5i H \-5 n < diam{P , . . . , P n }. 

Now, let L be any line hitting the disc -Br(o) and let / be the set of all 
indices i < n such that S% hits L. Applying the fact to the points Pi = L D 
and (5j = \Js/2^/el, i G 7, we get 

Let {ei,e2} be an orthonormal basis of M 2 and consider the two grids of 
lines Lj = e\ + (is/y/2)e\ and L 2 = e^ - + [isj\/2)e2, where i are integers 
with |i| < \/2R/s and ej~ are the coordinate axes with normal vectors ej, 
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j = 1,2. Note that each segment Si hits at least one of the two line grids. 
From the above fact, we get that 

/— n 
^E^<2 E \ 1 (LlnB R (o))<2(2R + nR 2 /( S /V2)) 

i=1 \i\<V2R/s 

< 2s~ 1 ttR 2 (-— + V2] < hs~ l irR 2 , 
\tt R J 

hence 

n 

YV^< 5V2s-V\R 2 < Ss^nR 2 . 

Lemma 4.3. Let A be a nonempty compact set in M. 2 , x,y two critical 
points of its distance function with distances < r < s from A and let a,b 
be points of A with \a — x\ = Aa(x) = r and \b — y\ = A J 4(y) = s. Then, 
\x — y\ > \/2r^s — r. Furthermore, if S (T) is the segment with endpoints x 
and (resp., y and then dist(5 , ,T) > ^/r~j2y r s~^? '. 

Proof. The interior of the ball B s (y) contains no point of A. Thus, 
int-B s (y) contains no diameter of B r {x) (since, otherwise, x would not be 
a critical point). It follows that \x — y\ > y/2r\/s — r (see [5], Lemma 4.3). 
Furthermore, since a G dB r (x)\vak,B s (y) and b G dB s (y) \ int B r (x), we have 
(6 — a) ■ (y — x) > 0. It follows that the orthogonal projections of the segments 
S and T onto the line passing through x and y have distance at least \y — 
x\/2 > y/r ' j2\j ' s — r. The same lower bound must then apply for the distance 
of the segments themselves. □ 

Corollary 4.1. Let iCR 2 be compact with D = diamA and letO < s <r. 
Write 

( S ,r)\C(A)=U/, 

i 

as a countable union of maximal disjoint open intervals with lengths £j. 
Then, 

i 

We now consider a Brownian motion. Given r > 0, we introduce the ran- 
dom variable 

rj r = inf{e > : r + e E C(S)}, 

where we use the convention inf = 00. By [12], Theorem 4.1, any r > does 
not belong to C(H) almost surely. Hence, we get 7] r > almost surely. For 
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any curve 7 : [0, 1] — > Mr, define its supremum norm by ||7|| = max^rgu |7(i)|- 
When we write \\X\\ in the following, we implicitly restrict X to [0, 1]. 

Lemma 4.4. Given r > 0, there exists a 5 > such that for all s S (r — 
5, r) and all e < 5, we have 

r(vs<(^) 3 \\x\\^> 1 -r(i 1 r<\\xfe). 

Proof. For each trajectory of X, choose a critical point z£l 2 with 
Ae(z) = r + r) r and denote the random variable £ = \X(1) — z\ — (r + rj r ). 
Since X(l) does not lie on the boundary of B r+Vr (z) almost surely (see [12], 
Lemma 4.2 or [9], Lemma 3.3), we have £ > a.s. Thus, as r \ 0, it holds 
that 

{Th.<\\X\\ i e,Z>T}/ , {r h .<\\Xfe}. 
By choosing r > small enough, we have 

P(7 ?r <||X|| 4 e ,e>r)>|P(7 7r <||X|| 4 e). 
Furthermore, we can choose Ao > 1 such that 

P( sup \X(t)-X(l)\<r)>j. 

Since a Brownian motion has independent increments, we have 

p(7? r <||X|| 4 e,£>r, sup \X(t)-X(l)\<r) 
V l<t<A 7 

= P(r/ r < ||X|| 4 e,e>r)p( sup \X(t) - X(l)\ < r) 

V l<i<A ' 

>Ap (r?r< ||x|| 4 £ ). 

Note that if rj r < e, £ > r and sup 1<t<A() |AT(i) — AT(1)| < r, then r + r/ r 
remains a critical value of {AT(i):0 < t < Ao}. By the scaling invariance 
property, E* = A~ 1/2 {X(i) : < t < A} has the same distribution as S, with 
critical value A _1/2 (r + 7/ r ) for any 1 < A < Ao- Denoting by ||X||*,r?* the 
random variables corresponding to H* [thus, ||AT||* = A~ 1/2 maxo<t<A |X(i)|], 
we have 77* 1/2 < A _1 / 2 7/ r , ||AT||* > A~ 1/2 ||X||. Consequently, we have 

ip( ??r < ||X|| 4 e) < A P(r/r < || X ||4 e) 

<P(?? r < ||Al 4 e,£>T, sup |X(t) — X(l)| < t) 
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= Pfe-i/2 r <A- 1 / 2 ||Al 4 e,£>T, S up \X(t)-X(l)\<T 

v Kt<A 

<p(T,*_ 1/ar <A 3 / 2 (||X|r) 4 e,e>r, sup \X(t) - X(l)\ < r) 

= p( Vx - 1/ 2 r <X 3/2 \\X\\ 4 e^>T, sup \X(t)-X(l)\<r 
v Kt<A 

<P(^- 1/2r <A 3 / 2 ||X|| 4 £ ) 



— 1/2 

for all 1 < A < Ao- The assertion thus holds with o = r — A r. □ 

Proposition 4.1. For any r > 0, the mean number of connected com- 
ponents of K 2 \ S r is finite. 

Proof. Let N r be the number of bounded connected components of 
M 2 \ S r . For r > 0, consider 7? r . If r\ r = oo, then H r does not have any holes 
and hence N r = [otherwise, any hole would contain a maximal inscribed 
circle with its center at a critical point and with critical radius s G (r, oo)]. 
Since r ^ C(H) a.s., we get C(E)n [r,r + 7y T .) = a.s. By Lemma 4.1, we thus 
have 

7T (2||X||+2r) 2 
N r <^= — a.s., 

where we adopt the convention that — = 0. We shall show that E||X|| 2 r/ r 1 ^ 2 
is finite, which is equivalent to the assertion. Let a = We have 



poo poo 

E\\X\\ 2 rj- 1/2 = / P(\\X\\ 2 r]- 1/2 >t)dt= P(r? r < ||X|| 4 r 
Jo Jo 



\dt 



2 



oo 



-X E / 1 (r) s <<T 3 ||X|| 4 i- 2 )^ s ^> 

JO Jr—8 

by Lemma 4.4 and Fubini's theorem, since 



i<JPfo. < \\X\\ 4 e) = \f P(r? r < \\X\\ A e) ds 

Z Z Jr—6 

< f P(t] s <a 3 \\Xfe)ds 

Jr-5 



E / ^{ns«J 3 \\x\\^)ds 



for any e > 0. Fix a trajectory H of X. Let {ij} be the family of maximal open 
intervals of the complement of C(E) which intersect the interval (r — 5, r) 
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and let £j be the length of We assume, without loss of generality, that 
the sequence {e.;} is nonincr easing. We then have 

/ 1 (r, s <a3\\x\\it-2) ds <J2 min{cr 3 1| A"|| 4 t- 2 , ej 

Jr-8 i 

and hence, by Corollary 4.1, 

E||X||V 1/2 < -E J ^mm{a z \\Xfr 2 ,Ei}dt 

= -EV / Eidt+ , r 2 dt 

= ^^(a 3 / 2 ||X|| 2 ^ + cx- 3 / 2 ||X||- 2 ^) 

j 

< ?E(a 3 / 2 ||X|| 2 + a- 3 / 2 ||X||- 2 )8(r - 5)- 3 / 2 ^||X|| 2 
o 

(a 3 / 2 E||X|| 4 + ^ 3 / 2 ). 



5(r - <5) 3 / 2 

Since all moments of ||S|| are finite, the proposition is proved. □ 

Corollary 4.2. The mean Euler number of the Wiener sausage H r in 
M 2 is finite for any r > 0. 

5. Convergence of parallel neighborhoods. Let A C M. d be a nonempty 
and compact subset of M d . For r > 0, denote by Va(?*) the volume Vrf(j4 r ) of 
the parallel neighborhood to A. The following result was proven by Stacho 
[17], Theorem 3. 

Theorem 5.1. Let {A n } be a sequence of nonempty compact subsets of 
M. d converging with respect to the Hausdorff metric to a nonempty compact 
set A. Then, 

(i) VA n {r) — > Va{t), n — ► oo, for all r > 0; 

(ii) H^idAlf) H^idAr), n — > oo , for all r > where (Va)'(V) exists. 

Remark 5.1. In fact, assertion (i) is included in the proof of [17], The- 
orem 3. It follows from the fact that for any e > 0, A r - £ C (A n ) r C A r+£ 
for sufficiently large n and from the continuity of the volume function Va- 
Assertion (ii) is formulated in a different form in [17], Theorem 3, but (ii) is 
exactly what is proved there. 
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We shall now show that the curvature measures (and, hence, also the 
intrinsic volumes) of the parallel neighborhood also converge. 

Theorem 5.2. Let {A n } be a sequence of nonempty compact subsets of 
]R rf converging with respect to the Hausdorff metric to a nonempty compact 
set A and let r G (0, oo) \ C(A). Then, r ^ C(A n ) for sufficiently large n and 
the curvature measures Ci(A^;-) converge weakly to Cj(A r ;-), n^oo, for 
i = 0,...,d-l. 



The proof will be based on the following result. 



Theorem 5.3 ([3], Theorem 5.9). Let D n ,D be compact sets with posi- 
tive reach such that dH{D n , D) — > and inf n reach D n > 0. Then, Ci{D n ; ■) — > 
Ci(D; •) weakly, n — > oo, i = 0, . . . ,d — 1. 

To prove Theorem 5.2, it is sufficient to show that 

lim inf reach C" > 0, 

n— >oo 

where C™ = M d \ A™. We shall start with some continuity results. 

We say that a function F from a metric space X to the space KL' of 
nonempty compact subsets of IR^ (equipped with the Hausdorff metric) is 
upper semicontinuous if, whenever x n — > x in X, a n — > a in IR'' and a n G 
F(x n ), n G N, we have a G F(x) (see [11]). Recall, also, that a real function 
/ : X — > M. is upper semicontinuous if lim inf f(y) < f{x) for all x G X. 

Lemma 5.1. (i) The function S^(x) is upper semicontinuous on 
R d for any compact A C R rf . 

(ii) The function A \— > S^(x) is upper semicontinuous on JC' for any x G 
R d . 



Proof. To verify (i), let x n — > x, a n G T,A(x n ), a n — > a. We shall show 
that a G S^(x). If not, there would be another point b G A with |6 — x\ < 
\a — x\. Let n be sufficiently large that 

max{|a n — a|, \x n — x\} < e = -(|a — x\ — \b — x\). 

3 

Then, by the triangle inequality, 

\x n — b\ < \x n — x\ + |6 — x\ < e + \a — x\ — 3e 
< \a — x\ — \a — a n \ — e < \a n — x\ — e 

^ I &n | ^ | — | \ j 

which means that a n is not the closest point of A to x n , a contradiction. 
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In order to show (ii), let A n ,A be compact sets, dn{A n , A) — ► 0, a n G 
S^n(a;) ) a n — > a. We have to show that a G ^^(x). If not, there would be a 
point 6 G A with |6 — a;| < \a — x\. By the definition of Hausdorff distance, 
there is a point b n G A n with |6 n — b\ < e for sufficiently large n, where 
e = |(|a — x\ — \b — x\). We also have \a n — a\ < e for large n. We then easily 
obtain, by the triangle inequality, that \x — b n \ <\x — a n \, which contradicts 
the fact that a n is the closest point of A n to x. □ 

As in [12], we introduce the function 

J A '• x i — ► min{|a — x\ : a G conv x G M. d . 

Clearly, x is a regular point of if and only if Ja(x) > 0. 

Lemma 5.2. (i) The function x i— > Ja(^) upper semicontinuous on R rf 
for any compact iCR 1 *. 

(ii) TTje function A i— ► Ja(^) *s upper semicontinuous on K! for any x G 

Proof. This follows directly from Lemma 5.1 and from the fact that 
Ja^e) depends continuously on T,a(x). □ 

Lemma 5.3. IfA n ,A are nonempty compact subsets o/M rf , d,H{A n ,A) —> 
as oo and r ^ C(A) ; then 

liminf inf{ J^fx) : x G dA™} > 0. 

n— >oo 

Proof. We have inf{ Ja{%) '■ x G dA r } > by [12], Lemma 3.3. Further, 
since J a is upper semicontinuous and dA is compact, we have 

rj = inf{ Ja(x) : r — e < Aa(x) < r + e} > 

for some sufficiently small e > 0. By the definition of the Hausdorff metric, 

— £ C A^ C A r j rE 

for large n. Now, for any x with r — e < Aa(x) < r + e, we have Ja™(^) > f?/2 
for sufficiently large n, by Lemma 5.2. Define 

K n = {x:r — e < Aa(x) <r + e, Ja™{x) < n/2 for m > n}. 

By the upper semicontinuity of Ja™, the sets K n are compact and since 
K n \ 0, we have if no = for some no, which completes the proof. □ 

Lemma 5.4. If A C R" is a nonempty and compact subset of M. d and 
r > 0, then 

re&chW d \A r > inf{ J A {x) : x G dA r }. 
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Proof. If z G dA r , then, for any a G convE^t^), B\ z _ a \(a) ^ with 

|^ — a\ > 7] = inf{ Ja(x) : x G cL4 r }. 

We can thus "roll" a ball of radius rj from outside along the boundary of 
R rf \ A r , which proves that reach M d \A r >w (see [3], Theorem 4.18). □ 

Proof of Theorem 5.2. Follows directly from Lemmas 5.3 and 5.4. 

□ 

6. Approximations of the Wiener sausage by polyconvex sets. In this 
section, we use a piecewise linear approximation H n of the path S of a Brow- 
nian motion X to construct the almost sure approximation of the Wiener 
sausage H r = H © B r (o) by random polyconvex sets 3™ = E n © B r (o). Under 
some assumptions, the convergence of the corresponding curvature measures 
and their first moments is shown in Corollary 6.1 and Proposition 6.1. In 
Sections 6.2 and 6.3, examples of such approximations are considered in 
detail. 

6.1. General convergence results. Let H n be an arbitrary random piece- 
wise linear curve with vertices X n (ti), i = 1, . . . , k n , approximating the path 
H of the Brownian motion in the sense of the Hausdorff distance: dn (E n , H) — > 
as n — > oo, almost surely. For each ii£N, consider its parallel neighbor- 
hood = E n © B r (o), r > 0. It holds that H™ G TZ a.s. for all n G N and 
r > 0. We call the sequence {H"} n6 pj of polyconvex random closed sets an 
( almost sure ) polyconvex approximation of the Wiener sausage E r . 

Since S" is polyconvex, its curvature measures Cj(H"; •) and intrinsic vol- 
umes Vi(H"), i = 0,...,d, are well defined for all r > 0. In the following, 
we show that the curvature measures of H" approximate the corresponding 
curvature measures of the Wiener sausage for growing n. 

The next result follows immediately from Theorems 5.1, 5.2 and [12], 
Corollary 4.1 and Lemma 4.3. 

Corollary 6.1. Let d//(S n ,H) ^0 as n— > oo, almost surely. We then 
have the following: 

(i) V=«(r) — > V=(r) as n — > oo almost surely for any r > 0; 

(ii) 7Y d_1 (<9H") — > 7l d ~ l (9H r ) as n — > oo almost surely for any r > if d< 
3, and for almost all r > if d>4; 

(iii) Cj(H"; •) — > Cj(H r ; ) weakly as oo almost surely for any r > and 
i = 0,...,d-l ifd<3. 

Corollary 6.1 combined with the dominated convergence theorem yields 
the following proposition. 
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Proposition 6.1. Let d h (H n , S) — >-0 as oo, almost surely. We then 
have the following: 

(i) EVs™(r) — > EV=(r) as n — > oo for any r > 0; 

(ii) WH d ~ 1 (d'E^) — ► KH d ~ 1 (dE r ) as n — > oo /or any r > if d<3, and for 
almost all r > if d>4. 

Proof. It suffices to find the uniform integrable random upper bounds 
for |V*(E?)|, i = d,d-l. 

(i) First, note that 

iii^i -imn<^(H™,H). 

Since the last expression tends to as n — > 0, we have 

V d (E?) < V d (B r+m (o)) =uj d (r + \\X n \\) d < u d (r + \\X\\ + l) d 

for sufficiently large n (and all u). It is known that all moments of the 
maximum of the Bessel process are finite; see, for example, [12], 

proof of Theorem 2.1. Consequently, we have 

EV d (Z?) < uj d E(r + \\X\\ + l) d <oo 

for all sufficiently large n. 

(ii) Since <9E™ C A^l({r}) holds for all r > 0, we use the co-area formula 
from [4] to obtain 

H d - 1 (Z?)<H d - 1 (A E U{r})) = Vi i (r) 

for all r > such that Vg (r) exists. These are all r ^ C(3 n ). By Theo- 
rem 5.2, it holds that r ^ C(E„) for all sufficiently large n if r ^ C(E). 
It is known (see [12]) that r ^ C(H) almost surely for any r > if d < 3, 
and for almost all r > if d > 4. By [12], Lemma 4.6, we obtain 

H d -^<^(r)<~^(i? + r) d , 

where R = \\X n \\ < \\X\\ + e for all sufficiently large n. Since the upper 
bound 

^ d (\\X\\+r + e) d 
r 

is integrable, the dominated convergence theorem yields EW« , - 1 (5H?) 
-> EH d_1 (9E; r ) as n -» oo. □ 

The following statement follows directly from Corollary 6.1(iii) and the 
dominated convergence theorem. 
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Remark 6.1. Let d#(H n , 3) — > as n — > oo, almost surely, d < 3. If the 
sequence of random variables {^(H™)} ng N is uniformly integrable for i = 
0, . . . , d- 2, then EVi(E?) -> E^(H r ) as n -> oo for i = 0, . . . , d- 2 and r > 0. 
We conjecture that in the two-dimensional case, the sequence {Vo(H")} ne pj 
is uniformly integrable and that convergence for the mean Euler number 
holds. 



6.2. An example. Consider a regular lattice of time moments ti = i/k n , 
i = 1, . . . , k n , k n — > oo as n — > oo. Let X n be the piecewise linear approxi- 
mation of X passing through X n {ti) = X{ti) for alH = 1, . . . , k n and n € N. 
It means that the vertices of the piecewise linear approximation X n of the 
Brownian motion X lie on the path H n of the Brownian motion. In the 
following, we show that such an approximation procedure preserves the con- 
vergence of the mean curvatures in the sense of the last section. 

Theorem 6.1. It holds that 

d H {Z n ,E) < \\X n -X||, n£N, 

almost surely. 

Proof. Using the piecewise linearity of H n , we get 
maxA Sn (y)= max max A Sn (X(t)) 

yer. «=l,-,fente[(l-l)/fen,V fe n] 

= max max min min \X n (s) — X(t)\ 

%=!,.. .,kn tG[(i-l)/k n ,i/k n ] i=l,...,k„ s£[(i-l)/k n ,i/k n ] 

< max max \X n (t) - X(t)\ = max \X n (t) - X(t)\. 

i=l,...,k n tel(i-i)/k n ,i/k„] te[o,i\ 

Using similar arguments, it can be shown that 

maxA H (y)< ||X n -X||. 

By the definition of Hausdorff distance, we have 

d H (E n ,E) =max(maxA s (j/),maxA=»(i/)} < \\X n - X\\. n 

Corollary 6.2. <i#(H n ,H)^0 as n^oo, almost surely. 



Proof. It is known that || W™ — Wi\\ — > as n — > oo almost surely for all 
i = 1, . . . , d. The assertion of the corollary follows from the above relation, 
Theorem 6.1 and the inequality 



\X n -X\ 



max 

te 



[0,1] \ 



J2(wm - wm 2 < y: \\w? - Wi\\ - o. 



□ 
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6.3. Haar-Schauder approximation. The piecewise linear approximations 
considered above can easily be simulated using Gaussian finite-dimensional 
distributions of a Brownian motion. Another constructive approach is pro- 
vided by the so-called Haar-Schauder approximation of a Brownian motion; 
see, for example, the book [6]. 

Let H k : [0, 1] — > R, k G N, be the so-called Haar function defined by the 
relations 



#l0) = l, a €[0,1], 



H 2 ™ + k{s) — < 



2m/2 

2 m /2 

10, 



s G 



fc-1 2Jfe-l 



2"i ' 2 m +i 
2fc- 1 fc 

2»n+l ' 2"i 



otherwise, 
The Schauder function is given by 



for k = 1, 2, . . . , 2 m and m = 0, 1,2, . 

S k {t)= f H k {s)ds, ken. 
Jo 

It is known that the Wiener process is equal in distribution to the Haar- 
Schauder series 

oo 

W(t) = Y / Y k S k (t), i€[0,l], 
fc=i 

where {l^i}neN is the sequence of i.i.d. iV(0, l)-distributed random variables; 
see [6], pages 56-59. This series converges (a.s.) absolutely and uniformly on 
[0,1]. Hence, the Wiener process can be approximated pathwise by partial 
sums 

2" 

Y,Y k S k (t), t€[0,l]. 
k=l 

We use this idea to approximate the coordinates Wi(t), i = 1, . . . ,d, of the 
Brownian motion X(t) a.s. by 



k=l 

where the sequences {Y ik } k& ^, i = I, . . . ,d, of i.i.d. standard normally dis- 
tributed random variables are independent. For 

X n (t) = (W 1 n (t),...,W2(t)), t€[0,l], 

consider its path H n = {X n (t) : < t < 1}. It is a piecewise linear curve with 
2 n + 1 nodes lying on S. More precisely, we have X n (0) = X(0) = o a.s. for 
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(a) Estimated Euler number 

Fig. 1. Mean (solid), minimal (dashed) 
and area, respectively, of 1000 simulated 
discretization steps. 
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n G N and 

X n (t) = (1 - a n (t))x(^Pj + a n (t)x(^) , t G 

a.s. for all n and k = 1, . . . , 2 n , where a n (t) = 2 n t — k + 1. Thus, 5 n can be 
regarded as an a.s. approximation of the path 3 by piecewise linear curves, 
in the sense of Section 6.3, with k n = 2 n . 

7. Numerical results. In this section, we present some numerical results 
for the estimated mean intrinsic volumes of the Wiener sausage in the two- 
dimensional case. The theoretical basis for the simulation which supports 
the outcome is provided by the convergence results of Proposition 6.1 and 
Corollary 6.2. 

7.1. Approach and estimation results. A piecewise linear approximation 
of the path of a Brownian motion in R 2 is considered and the path is dilated 
to produce a polyconvex approximation of the Wiener sausage. The approx- 
imation of the path is achieved by simulating the independent, normally 
distributed increments of the Wiener process; see Section 6.2. For this, we 
chose the time interval [0,T] = [0, 10 6 ] and k equidistant time points in this 
interval. Figure 1 shows the estimated Euler number and estimated area, 
respectively, of approximated Wiener sausages 3 r in the plane with radius 
r = 20 depending on the number k of such discretization points. For each k, 
one can see the mean (solid) as well as the minimal (dashed) and maximal 
(dotted) values of 1000 realizations. It is obvious that for small values k, the 
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approximation is not very good and the estimated value of the polyconvex 
approximation is therefore too small. On the other hand, both estimated 
values converge to a limit once k is sufficiently large. Since this convergence 
behavior is basically the same for the boundary length, we discretized the 
interval [0,T] in k = 100,001 time points for all simulations. The radius r of 
the Wiener sausage varies between 1 and 1900 pixels. For radii less than 1, 
there is no difference between the Brownian motion and the Wiener sausage 
in the simulation due to discretization. In any case, in Figures 2(a)-2(f), 
one can recognize a trend in the curves for small radii. On the other hand, 
the behavior of the curves for growing radii is clearly visible, so there is no 
need to simulate bigger radii. In Figures 2(a)-2(f), the solid line denotes 
the mean value, the dashed line the minimal value and the dotted line the 
maximal value of the considered estimator out of 1000 realizations for each 
radius. The intrinsic volumes were estimated with the algorithm described 
in [7] using Steiner's formula. Figure 2(a) shows the mean, minimal and 
maximal estimated area of 1000 realizations. For growing radii, the area 
grows quadratically and as the radius tends to zero, the area tends to zero 
as well. This can be seen in Figure 2(b), which shows a magnification of 
Figure 2(a) for small radii. Thus, we get an empirical illustration of the 
following asymptotic result of Le Gall [10]: 

ttT 

EV=(r)~- r , r^0. 

I log 7" | 

Figures 2(c) and 2(d) show the estimated mean boundary length. One 
can see that it grows linearly once all holes of E r disappear as r tends to 
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Fig. 2. Estimated mean intrinsic volumes 
depending on the radius. 
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(c) Mean (solid), minimal (dashed) and 
maximal (dotted) boundary length of 
1000 simulated Wiener sausages depend- 
ing on the radius 
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(e) Mean (solid), minimal (dashed) and 
maximal (dotted) Euler number of 1000 
simulated Wiener sausages with large 
radii (r > 120) 




(d) Magnification of Figure 2(c) for small 
radii 




(f) Mean (solid), minimal (dashed) and 
maximal (dotted) Euler number of 1000 
simulated Wiener sausages with small 
radii (r < 15) 



Fig. 2. Continued. 



infinity. It is shown in [12] that the expected boundary length equals the 
first derivative of the mean area of the Wiener sausage with respect to its 
radius. This is also backed by our simulations. For instance, the gradient 
of the tangent line to the graph of the mean area in Figure 2(b) grows 
rapidly as the radius becomes small. This coincides with the behavior of the 
graph of the mean boundary length in Figure 2(d) around the origin and 
the asymptotic formula 

ttT 

KH 1 (dE r )~-—z-, r^0, 
rlog r 

proved in [12]. 
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The mean Euler number of H r , that is, 1 minus the mean number of 
holes, tends to — oo and its slope tends to infinity as r — > since H r is 
connected and the number of holes grows unboundedly. This behavior is 
reflected in Figure 2(f). For growing radius, the mean Euler number tends 
to 1 as more and more holes get filled by the dilation, and equals 1 when all 
holes disappear; see Figure 2(e). 

Berezhkovskh, Makhnovskh and Suris [1] give an explicit formula for the 
expected area of the Wiener sausage in dimensions greater than 1. Its special 
case for d = 2 can be written as 

EV E (r) = 7rr 2 + — i * - ■ ^ dy, 



* Jo y*(JZ(y) + Y2(y)) 

where Jo and Yq are Bessel functions of order zero. Thus, EV=(r) and its 
derivative with respect to r — the expected boundary length WH l (d'B r ) — 
can be computed numerically, in addition to the simulations shown above. 
Our computations show that the difference between the simulated and the 
approximated values is less than 10% for radii greater than 20 (resp., 200) 
for the boundary length (resp., area). For smaller radii, the difference gets 
bigger (up to 39%) due to the high slope of the boundary length and area, 
respectively. 

7.2. Approximation formulae. The curves of the estimated mean intrin- 
sic volumes were approximated with the Curve Fitting Toolbox in MATLAB 
using a trust region algorithm. The mean estimated area [see Figures 2(a) 
and 2(b)] of the Wiener sausage can be approximated by the function 

EVs (r) « + 387768r ' 252344 + 1064.51r 1 ' 24465 

v ; | Iogr| + 265.265 

+ 0.235594r 2 ' 27655 

depending on the dilation radius r. The absolute value of the relative error 
between this function and the estimated values is always less than 2% and 
the mean relative error is less than 0.6%. 

For the boundary length in Figure 2(c), we fitted the function 

EHHdZr) « 5 -^^ + 93484.2r-°' 862668 

v ; r(logr) 2 + 2.66290 

+ 9038.097- ' 123680 + 15.64 93r a9019n . 

This is not exactly the derivative of the approximated mean area, but, nev- 
ertheless, the fit is very good. The relative errors are less than 1.9% for all 
radii and the mean relative error is less than 0.6%. In both cases, that is, for 
the estimated mean area and boundary length, their asymptotic behavior as 
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r — ► coincides (up to a constant factor) with the formula of Le Gall given 
in the previous section. 

A different class of functions was used to fit the Euler number. Here, a 
logarithmic term seems to be a good choice to approximate 2 — EVo(H r ). 
[We did not approximate EVo(H r ) directly because it takes values close to 
zero in a certain interval and therefore one cannot compute relative errors 
to evaluate the fit.] Since this logarithmic function converges to 1 slowly 
as r — ► oo, we additionally used the cumulative distribution function $ of 
the standard normal distribution to make it converge faster so that it is 
practically 1 for large radii. The resulting function fitted with MATLAB is 
given by 

¥vr) 0.0423017(1 - $((r-224.899)/50.2096)) 

EV {^r) « " l og (3. 88182 . l -6 r 1.88978 + L0 ) r 0.153452 ■ 

The relative error of this fit is less than 10.5% and the mean relative error 
is less than 1.5%. For growing radii, the fitted function converges to one 
at an appropriate speed, which leads to low approximation errors. For radii 
greater than 160, the relative error is less than 4% and for radii greater than 
560, the relative error is less than 1%. 

Acknowledgment. The authors are grateful to Volker Schmidt for fruitful 
discussions. 
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